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Abstract. A ID Dirac tight-binding model is considered and it is 
' shown that its nonrelativistic limit is the ID discrete Schrodinger model. 

For random Bernoulli potentials taking two values (without correla- 
tions), for typical realizations and for all values of the mass, it is shown 
that its spectrum is pure point, whereas the zero mass case presents 
dynamical derealization for specific values of the energy. The massive 
case presents dynamical localization (excluding some particular values 
of the energy). Finally, for general potentials the dynamical moments 
for distinct masses are compared, especially the massless and massive 
■ Bernoulli cases. 
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Besides the huge amount of mathematical works on spectral problems 
_C ■ related to the one-dimensional Dirac model |IJ in physics it has also 

been used in comparative studies of relativistic and nonrelativistic electron- 
localization phenomena in relativistic investigations of electrical con- 
duction in disordered systems in the construction of supertransparent 
models with supersymmetric structures and in relativistic tunnelling 
^ ' problems |2Uj . 



1. Introduction 



■ In this paper a discrete version of the ID Dirac model is discussed, which 

can be interpreted as a relativistic version of the well-known tight-binding 
Schrodinger Hamiltonian (with h = 1) 

(1) {Hlp) n = - — (AV>)n + V n l/> n = —(-1p n+ i - 1p n -i + 2lp n ) + V n 1p n . 
2m 2m 

The model was first reported in [11 and this work is its very expanded 
and mathematical detailed version. Consider a particle of mass m > in 
the one-dimensional lattice 7L under the real site potential V = (V n ). The 
proposed ID Dirac tight-binding operator is 

(2) K)(m,c) = TD (m,c) + V Id 2 = cB + mc 2 a 3 + V Id 2 , 
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with c > representing the speed of light, 

B , 



d 

a"3 the usual Pauli matrix, Id<i the 2x2 identity matrix and d the finite 
difference operator (a discrete counterpart of the first derivative) defined by 

(dl/)) n = 1p n+ l - Ipn- 

(d*ifi) n = V'n.-i — ipn i s the adjoint of d so that JDo(m, c) = cB + mc 2 az is a 
bounded self-adjoint operator acting on 1 2 {7L\ (D 2 ) and its square is 

„ , v2 / c 2 o!d* + m 2 c 4 

= V cfdtT + mV 

This equality is reminiscent of the relation between momentum p and en- 
ergy E in relativistic quantum mechanics [I], given by E 2 = cPp 2 + m 2 c 4 . 
Denoting by cr(A) the spectrum of a self-adjoint operator A, it is well known 
that cr(-A) = [0,4], and since d*d = dd* = -A, 



cr(IDo(m, c)) = — cy 4 + m 2 c 2 , —mc 2 U mc 2 , cy 4 + m 2 c 2 

In case the potential V is a bounded sequence, ID(m,c) is also a bounded 
self-adjoint operator acting on £ 2 (2;(D 2 ). 

It will be shown that the nonrelativistic limit of the resolvent of the dis- 
crete Dirac operator (J2J) is the resolvent of discrete Schrodinger operator (jTjl 
(when projected on a proper subspace; see Section EJ). This is an important 
support for such Dirac model. 

The study of quantum transport depends, of course, on the admitted defi- 
nitions. In the physics literature terms like "extended states" and "zero Lya- 
punov exponents" have been used to crudely designate quantum transport. 
For instance, in |^7j it was claimed that "extended states" were found in one- 
dimensional Schrodinger systems with off-diagonal randomness, but in [21] 
it was argued that although the localization length diverges the "transmis- 
sion coefficient" vanishes as the system size goes to infinity. Up to recently, 
in the mathematical literature pure point spectrum (sometimes with expo- 
nentially decaying eigenfunctions) was considered synonymous of absence of 
transport. Currently the transport has been probed via the time behavior 
of the moments of the position operator, and in this work this idea will be 
followed. See ahead for precise definitions and related comments. 

One of the goals of this paper is to study the phenomenon called dynamical 
localization (in the sense of time-boundedness of all moments of the position 
operator) for the Bernoulli-Dirac model, that is, the model (J2J) with the site 
potentials V n , n £ 7L, being independent identically distributed Bernoulli 
random variables taking the values ±V, V > 0. In this case it will be shown 
that almost surely the spectrum of ID(m, c) is pure point for all values of the 
mass, the massive case has dynamical localization (excluding some particular 
values of the energy for which a more careful analysis is needed) and the zero 
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mass case presents dynamical delocalization (that is, absence of localization) 
for specific values of the energy. 

The problem of dynamical localization has been intensively studied dur- 
ing last years, especially in the case of random discrete and continuous 
Schrodinger operators (in particular for the Bernoulli- Anderson model, that 
is, the Schrodinger model with Bernoulli potentials); see ^Ql EE1 EE! and 
references there in. What one usually proves is the so-called exponential 
localization (UIEIEHI, i-e., pure point spectrum and exponentially decaying 
eigenf unctions. On the other hand, it is also known that exponential lo- 
calization does not imply dynamical localization ^21 > it is usually needed a 
precise control of the decay of the eigenfunctions, called SULE jl21ll(ij . that 
can be obtained through the method of multiscale analysis, a technique set 
out by Frohlich and Spencer |14l I15j . 

One motivation for studying dynamical localization for the Bernoulli- 
Dirac operator comes from the random dimer model |131llUj . i.e., the Bernoulli- 
Anderson model with the site energies V n assigned for pairs of lattices: 
V2n = V2 n +i = ±V for all n. This model almost surely presents pure 
point spectrum for all values of V 7^ [TJ3- It was also numerically found 
in ^3] and rigorously shown in |l(Jl I18j the existence of critical energies (in 
the sense of ^Hh see ahead) at which the Lyapunov exponent vanishes; dy- 
namical localization was obtained in 10 only after projecting onto closed 
energy intervals not containing such critical energies. Despite the similar- 
ity between the transfer matrices of the two models, it is not immediate 
the adaptation of the localization (delocalization) results to the Bernoulli- 
Dirac model and each step needs to be verified; here, many points will not 
be detailed when they follow exactly the same lines of their Schrodinger 
counterpart. 

With respect to nontrivial quantum transport, probed via dynamical de- 
localization (unbounded moments of the position operator) , it was found in 
random polymer models JI2, and in random palindrome models (both 
including the important random dimer model), due to existence of critical 
energies ^Hl- Recently, for ID discrete Schrodinger operators, Damanik, 
Siito and Tcheremchantsev (Hj have developed a general method which al- 
lows one to derive quantum dynamical lower bounds from upper bounds on 
the growth of norms of transfer matrices, and they applied this method to 
some substitution, Sturmian and prime models, among others. Damanik, 
Lenz and Stolz [Hj presented an extension of this method to ID contin- 
uous Schrodinger operators, with application to the continuous Bernoulli- 
Anderson model. Another method to obtain quantum transport from upper 
bounds on transfer matrix was lately developed by Germinet, Kiselev and 
Tcheremchantsev |17j . with application to Schrodinger operators with ran- 
dom decaying potentials, providing new examples of Schrodinger operator 
with point spectrum and nontrivial quantum transport. 

In the zero mass case, the one-dimensional Bernoulli-Dirac model pre- 
sented here has pure point spectrum and nontrivial quantum transport for 
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potentials with no correlations nor decaying properties (see Section EJ). This 
phenomenon does not take place in the corresponding Schrodinger tight- 
binding model |16j . and this also motivates the interest to better understand- 
ing the Dirac case. Presumably, this tight-binding model is the simplest one 
presenting such phenomenon. 

Now the localization results for the Bernoulli-Dirac model will be briefly 
summarized. By using as the main tool a particular form of Fursten- 
berg Theorem (Lemma ahead) , it is shown (see Theorems El and 0} 
that the Lyapunov exponent ^ m (E) is strictly positive for the energies 
E £ a(JD(m,c)), except for: E = ±V, V 6 (0 ,c], V ^ c/y/2, in the 
case m = 0; and if m > for (Ey = 0, V = cy/2 + m 2 c 2 ) and the four 
energy-potential pairs (E v = ±cV2 + m 2 c 2 ±c/\/2, V = c/y/ty. 

For all energies E for which ^y m (E) > 0, a initial estimate for localization 
(LemmaK and the Wegner's estimate (LemmaE)) will be checked; by adapt- 
ing the method multiscale analysis |28 | I15 [ ITT)] to this model, it will be shown 
(see Theorems 12 and 13} that for typical realizations the spectrum of ID(m, c) 
is pure point and the corresponding eigenfunctions are semi-uniformly ex- 
ponentially localized (SULE) |12[I16|. This and the results of 16 a (properly 
adapted to ID(m, c)) imply dynamical localization. 

In the massless (m = 0) case, the values E = ±V with V £ (0, c], V ^ 
c/V2, are critical energies for the operator ID(0, c) and this implies (almost 
surely) upper boundedness for the transfer matrices in the vicinity of these 
energies (see Lemmas El and EJ) . By adapting the ideas of ^5] (see also jH] ) 
to ID(0,c) it will follow (see Theorem EJ) that for an initial spinor \F well 
localized in space, there is < C q < oo such that 

r±e-*' T M®(0 t t)dt>C q T*- 1 ' i , 
Jo 1 

for almost all realization of the potential (or exponent q — 1 instead of q — 1/2 
for every realization), where (X is the usual position operator) 

i.e., there is nontrivial quantum transport despite the absence of a continu- 
ous component in the spectrum of ID(0, c). 

In the case of the set of pairs (E v = ±cV2 c 2 ±c/V2, V = c/y/2) and 
(E v = 0, V = cV2 

_l_ m 2 c 2j ^ j g s i lown ( see Theorem |1J) that the Lyapunov 
exponent 7 m vanishes, but it was not possible to give an answer about 
dynamical localization for them. Nevertheless, for these cases there is a 
general dynamical upper bound (in fact valid for all potentials V) established 
in Theorem El 

For distinct masses m,m' > 0, but m close to m' , it is expected that the 
moments M^\m,t) follow closely the moments M^\m',t) (both with the 
same potential), at least for a small period of time. The final result to be 
reported is an inequality confirming such expectative; by making using of 
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Duhamel's formula, it will be shown (see Theorem 0) that, for the initial 
state ^ with only one nonzero component, there exists K q > so that, for 
all t > 0, 

< K q \m - m'\c 2 t q+2 . 

In particular, for the Bernoulli-Dirac model this relation with m' = gives 
quantitatively an estimate of how, for small times, the dynamics of the lo- 
calized regime follows the delocalized one (see also Corollary ^ in Section 0). 

This paper is organized as follows: In Section |2] the nonrelativistic limit 
for the discrete Dirac model (J2J) is discussed. In Section El the results about 
spectral properties of such model, dynamical localization (derealization) 
and a dynamical upper bound for moments are presented, whose proofs ap- 
pear in Section [5J In Section ^ some tools used in those proofs are collected. 
Finally, in the Section the dynamical moments with different masses are 
compared; in particular the dynamics of the massless and massive Dirac- 
Bernoulli cases. 



2. Nonrelativistic Limit 

In this section consider ID(m, c) with m > fixed and c as a parameter. 
For simplicity, ID(c) will denote ID(m, c), which is supposed to be self-adjoint 
with (real) potential V . 

The nonrelativistic limit means c going to infinity, and since the rest en- 
ergy mc 2 is a purely relativistic quantity, (as usual) it must be subtracted 
before taking this limit. The norm convergence of the resolvent operators 
(lD(c) — mc 2 — z) 1 , for z G (D\IR will be considered. A below is the pro- 
jector onto the subspace of "positive energies," and so Ai^ corresponds 
to the Schrodinger operator (0). It is interesting to compare the approach 
presented here with the one in 

Theorem 1. If z € <D\]R, then 

lim (ID(c) - mc 2 - z)~ l = A (H^ - z)' 1 , 

c— >oo 

1 B 2 ~ 

where A = — (Jcfe + 03) and Hqo = h VA, and the limit is in the norm 

2 2m 
of bounded operators. 

Lemma 1. If z € (D\IR, then 

(3) (10(c) - mc 2 - zy 1 = (A + |±£) S(c) (id + V^S(c)} ' , 



M^(m,t)-M^(m',t) 
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where Id is the identity operator and 
(4) 

' Id 



S(c) = ^Hoo - z 
Proof. Note that 



2mc 2 



2mc 2 



(]D (c) + mc 2 + z) (IDo(c) - mc 2 - z) = c 2 B 2 - 2mc 2 z - z 2 . 



Hence 



(5) (T> (c) - mc 2 - z) 



-i ID (c) +mc 2 + z ( B 



2mc 2 

A + c ® + - ) s 

2mc 2 ' 



2m 



2mc 2 



( B 2 z 2 V 1 
with So = ( z — 2 j ■ O n the other hand, by using the operator 



relation 



\2m 



(A + B)- 1 = (Id - A- l B)~ 1 A~ 
B 2 z 2 

with A = z and B = VA, one obtains 

2m 2m 



S{c) = S (ld + VAS { 



(6) 



Therefore, by © and © it is found that 



(K(c)-mc 2 -z) 1 = (E>o (c) - mc 2 - z) 1 (id + V (JD (c) - mc 2 - z 



-i 



cB + z\ ( ~ ~cB + z n 

A + - — - I S [Id + VAS + V— — 5-5 



A + 



2mc 2 
cB + z 



2mc 2 



_ cE 4- z 
, S{c)(ld+V-- T S{c) 
Zmc^ I \ 2mc z 



-i 



□ 



Proof. (Theorem ^) Since (-ffoo — z) is bounded for z £ (D\IR and 

,.2 



2mc 2 

for c sufficiently large, one can expand 



(7) 



OO / 

S(c) = E ( 



n=0 



\2mc 2 



(H 



< 1 



(H c 



where the sum is convergent in the operator norm. 
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< 1 



2mc 2 



and so 
(8) 

Replacing Q and ( 
with 



(^+T( C )r i = ^(-r( C )r 



n=0 



into © one obtains the expansion 

" Rn(z) 



(m(c) 



mc — 2 



r'=E 



n=0 



i?o(z) =A(fT 0O -z)- 1 , 
K{H 00 -z)- l ^- + ^-{H 00 



A, 



and so on, and the sum is convergent in the operator norm. The result then 
follows. □ 



3. Localization Results 
Consider the family of Dirac operators 

(9) JDUm,c)=( m ^ ^a^+VM, coen = {-V,V}*, 

on ^ 2 (2Z;C 2 ), where V w (n), n G ZZ, are i.i.d. Bernoulli random variables 
taking the values ±V, V > 0, with common (nontrivial) probability measure 
fi and product measure P = FlnGZ A 4 (^>( n )) • Let ^Jm ^ e * ne s P ec tral 
projector of D w (m, c) onto the interval JcK. 

Denote by the elements of the canonical position basis of £ 2 (2;(D 2 ), 

for which all entries are f ^ ^ except at the nth entry, which is given by 
J J and ^ ^ J for the superscript indices + and — , respectively. If 
\P = I ^_ J is a solution of the eigenvalue equation 

then it is simple to check that 

V> + (n + 1) ^ _ ^(n)/^ ( V> + 0) 



ip (n) y m y ?/> (n — 1 
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with 



(l 



Tl{E) 



+ 



2 4 

m c 



mc 



E + V 



\ 



J 



The transfer matrix from site k to site n is introduced by 

T^(E;n,k) = T^ n - 1 \E)T^ n - 2 \E)...T^ k \E) , n>k 

and Tm(E; n, n) = ld.2- For q > 0, let jXl 9 be the moment of order q of the 
position operator on ^ 2 (2Z; (D 2 ), i.e., 

\X\^ = Y,\n\ q m^)5+ + (5-^)5-). 



Definition 1. The operator JD^fjn, c) is dynamically localized on a spectral 
interval I if for all q > and for all exponentially decaying initial state 
^ £ £ 2 (7Z;<E 2 ) 

supM$\ u (m,t) :=sup (Pf m e- iJD ^ m ' c)t ^, \X\ q Pf m e- m - {m ' c)t ^} < oo 
t ' ' * 

P almost surely (P-a.s.). Otherwise TDu(m,c) is dynamically delocalized 

on I. If I = a (ID vim, c)), i/ien ^(m, t) will be denoted by M^(m, t). 



It is important to notice that although the Dirac operator acts on spinors, 
its eigenvalue equation, in the transfer matrix form, looks exactly like the 
equation for a one-dimensional Schrodinger operator acting on scalar valued 
functions, with the transfer matrix being in SX(2,IR). Hence the methods 
used in studies of the usual one-dimensional Anderson model, as Fursten- 
berg's Theorem, can be applied for this Dirac model; see Sections 0] and 03 

The localization results are gathered in the following set of theorems. 



Theorem 2. Let (D u (m, c))^ e n be as in 0) andV £ (0,c],V / c/y/2. 
Then, P almost surely, the Lyapunov exponent 

lm {E) = lim ^ln\\T%(E;n,l)\\ 

n—>oo \ fi\ 

exists, is independent of uj, and 

(i) (i.l) -f m (E^±V)>0 form>0, 

(1.2) lm {E = ±V) >0form>0, 

(1.3) lo (E = ±V) = 0. 

(ii) Let m > 0; then P-a.s. a(W) u] {m, c)) is pure point. 

(iii) (iii.l) Let m > 0. Then P-a.s. the operator ID w (m,c) is dynamically 

localized on its spectrum. 
(hi. 2) For any closed interval I C o"(ID w (0, c)), with ±V /, the 
operator H\,(0,c) is dynamically localized on I. 
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Theorem 3. Let (ID w (m, c)) wgf2 be as in 0), m > and V > c, V / 
c\j2 + m 2 c 2 . Then, P-a.s. the spectrum of 1)^(771,0) is pure point and this 
operator is dynamically localized on its spectrum. 

Theorem 4. Let (ID^rra, c)) UJg Q be as in 0J, m > and V = c/y/2 (re- 
spectively V = cy2 + m 2 c 2 ). Then the same conclusions of Theorem 
(resp. Theorem^ hold except at the four possibilities of energies Ey = 
±c\/2 

_l_ vp 2 ( p,J slC j (resp. Ey = 0). (The point is that Ey = ±c\/2 + m 2 c 2 ± 
c/\/2 (resp. £\/ = Oj are energies so that <y m (Ey) = 0.) 

For the next result it is convenient to use the average dynamical moments 

(10) 4l(m,T):=J Q ±e-^ T Mi%(m,t)dt, 

defined for m > and T > 0. The main reason for working with this kind 
of Laplace transform average is relation Q14J1 ahead. 

Theorem 5 (massless case). Let (10^(0, c)) 4JgfJ 6e as in (0) and V G (0, c], 
V 7^ c/y/2. Then, for q > and ^ ura'ia onZy one nonzero component, there 
exists < C q (uj) < oo such that, for T > 0, 

(i) A^(0,r)>C,HT9-V2 p_ a . s ., 

(ii) 4 9 iL(0> T ) > C q {w)T*-l for every u, 

i.e., ID^(0, c) is not dynamically localized on its spectrum. 

The following theorem establishes very general upper bounds for the dy- 
namical moments of the position operator; notice that it holds for any po- 
tential sequence V and is not restricted to the Bernoulli case. 

Theorem 6. Let ID(m, c) be as in m > and ^ with only one nonzero 
component (so in the domain of \X\ q for all q > 0). Then for any q G IN 
there exists < K q (V ,m, c) < oo such that 

M^\m,t) < K q (V,m,c) t q , t > 1. 

Remark. It is possible to adjust the constant K q so that the above up- 
per bound holds for t > e for any given e > 0, instead of just t > 1. 
Since M^\m,t) > M$'\m,t) for q > q', it is evident that M^\m,t) < 
Kr q i (V, m, c) i^l for real q. 
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4. Tools 

In this section some tools and notations that will be used in the proofs of 
the results presented in Section |31 are collected. For studying the positivity 
of the Lyapunov exponent 7 m , m > 0, the following particular form of 
Furstenberg Theorem [5] will be used: 



Lemma 2. Let Q m (E) be the smallest closed subgroup of SL(2, 1R) generated 
by the matrices T^(E) and T^ V (E). Then j m (E) > if 

• Q m (E) is not compact, and 

• there is no probability measure on P(IR 2 ) (the set of all the directions 
ofTR 2 ) that is invariant under the action ofQ m (E), which is equiv- 
alent to the statement: the orbit Q m (E) ■ x := {T ■ x, T G G m (E)} of 
each direction x G P(IR 2 ) contains at least three elements. 



If L > 0, n£2, consider the finite subset of 7L 

L 



Ai(n) = | A; G 7L : \k - n\ < ^ 

with boundary 

dA L (n) = {(k, k') : k G A L (n), k' A L (n), \k - k'\ = 1}. 

Denote by ~\D^ L (m, c) the operator ID^m, c) restricted to £ 2 (A^(n); (D 2 ) 
with zero boundary conditions outside A^(n). 

The matrix elements of an operator O on £ 2 (2; (D 2 ) are given by 

V (6-,05+) {S-,OSj; 
with "norm" 

\\O nk \\ 2 = m,05+)\ 2 + K5+, 05^)\ 2 + \(8-,05+)\ 2 + \(5-,05^}\ 2 . 

Now two important results required for the multiscale analysis are de- 
scribed. The first one is the Wegner's estimate, adapted from |B] to the 
discrete Dirac operator (details will be omitted, since they are long and 
very similar to the Schrodinger case): 



Lemma 3. Let ID w (m,c) be as in J^) and I a compact energy interval. 
For any 8 G (0,1) and r > there exist Lq = Lq(I, 9,t, m) > and 
a = a(1, 6, r, m) > such that 

P jo; : dist (£,a(lD^( )(m,c))) < e^j < <r aL * 

for all E £ I and L > Lq. 
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The second result is the initial estimate for localization, adapted from |2H] 
(details omitted): 



Lemma 4. Let D w (m, c) be as in e > and 9 £ (0, 1). For each Eq G IR, 

Eq G" a(jD^ L ^ (m, c)) u>ii/i 7 m (£b) > i/iere ezisi Lo = Lq(Eq, e,8,m) > 
and r = v(Eq, e, m) > such that 

-i 



lo : 



ID A i (0) i 



m, c 



OA- 



< e -(lm(Eo)-e)L/2 y k G 3A L (0) I > 1-e 



/or a// L > Lq. 



In order to obtain dynamical localization from the multiscale analysis, 
the following properties of JD u (m, c) are useful: 

(PI) With respect to the spectral measure of ID^m, c), almost every en- 
ergy is a generalized eigenvalue, i.e., with polynomially bounded eigenvector 
(see H3HU). 

(P2) If E ^ a(JDt Un \m,c)) and * G ^ 2 (2Z; (D 2 ) so that D w (m,c)* = 
then 

¥(») = -(D^(»)(m, C )-£)^[ J ) - 1) 

-(ip^(»W)-^(° S)*(fe + i), 

with {(h,h - 1), (| 2> l 2 + 1)} = 0A L (n). 



Property (P2) follows after defining the boundary operator T\ L ( n \ by its 
matrix elements 




if j - 1 = k, j G A L (n), k <£ A L (n); 
if j + 1 = k, j G A L (n), k $ A L (n); 
otherwise, 



noting that Z 2 (2Z; (D 2 ) = Z 2 (A L (n); (D 2 ) Z 2 (Z \ A L (n); (D 2 ) and 
D w (m,c) = ^W(m )C )+D^)K C )-^ AlW . 



In the zero mass case (m = 0) the operators IL\,(0, c), a; G presents 
critical energies E v = ±V for V G (0, c] , V / c/>/2, as defined in [TH], 

since 

either T Q V (V) = Id 2 and T ~ V (V) is elliptic (that is, (trace T ~ v (V)\ < 2) or 



12 CESAR R. DE OLIVEIRA AND ROBERTO A. PRADO 

Tq V (— V) = Id,2 and Tq (—V) is elliptic. Thus there exists a real invertible 
matrix Q such that 

u y sin(jjj-) cos(r]±) J 

Since the eigenvalues of this matrix are e %ri± and e~ ir)± , for both of the above 
cases one has n+ — 77 _ ^ kir, k G 7L (a condition required in ^Hl)- By using 
the modified Priifer variables, phase shifts, oscillatory sums, large deviation 
estimates as in |18| . one obtains the following result 



Lemma 5 (massless case). Let X > be arbitrary. Then there are b > 
and C < oo such that for every N G IN, there exists a set f2Ar(A) C with 
P(n N (X)) < Ce~ bNX and 

\\Tg(E;n,k)\\ < C 

for all uj G n\n N (X), < k < n < N and E G [Ey - iV^ 1 / 2 ,^ + 
jY-A-l/2]_ 

On the other hand, since \\Q T^ v (Ey) = 1, expanding T^ v (Ey + e) 

into powers of e one obtains 

\\Q T± v (E v + e) Q^W < l + a\e\ 

for |e| < <5, 0<a<oo, and one deduces the following 



Lemma 6 (massless case). For 5 > there exists C < oo such that for all 
n,k eZ and E G [Ey - 5,E V + 5], 

\\Ttf(E;n,k)\\ < Ce csln - kl . 



An inductive argument shows that, for £ G (D and m > 0, 
(11) 

n-l 

T£(E + (; n, k) = T%(E; n,k)~Cj2 T ^ E + ^ V + 1) S C^ l ) T ^ l > k )' 

l=k 

where 

Now, for z G (D\IR and m > 0, introduce the two-components Green's 
function 



Ct, J-.:n) \ _ I iS+,(T> UJ (m,c)-z) 1 ^ 
J~, (E) w (m,c) - z) -1 ^ 
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so that 

(■w-».c)-»>(^$ ) = #<»>■ 

By using transfer matrices, one obtains for n < 0, 
(12) 



and for n > 1, 

G- (*;n-l) i mi ' J V G- w (z;0) 



(13) 



For z = E + i/T (T > 0) and m > 0, it is also valid the following identity 
(adapted from Lemma 3.2 in |19j): 

(14) Af (m,T) = ^-Y,\ n \ q f ^ M;n)\ 2 + \G-,„(z;n)\ 2 )dE. 



5. Localization Proofs 
In this section the proofs of Theorems [21 - ED are presented. 

Proof. ( Theorems 121 and 

The strategy of the proof is based on reference jTU], where the random 
dimer Schrodinger operator was studied. Since for the discrete Dirac op- 
erator there are the particular role played by the mass and some different 
possibilities for the transfer matrices, a rather detailed proof will be pre- 
sented. The idea is to show that given e > 0, / C a (JD w (m,c)) a com- 
pact energy interval not containing the excluded V values, then for each 
< 7 < 7m(-0 := inf{7m(-E') ■ E £ 1} there exist a constant C(uj,€, 7) and, 

f <pf\ 

for each eigenfunction ra,- w = ir \ with energy E~ w £ I, a "center" 

V <p,, u J 

Zj jU1 6 2Z, such that 

(15) ||^>(n)|| < C(a;,e,7)e 7l ^ |e e- 7|ri -^ 1 , Vra e Z. 

If \E r decays exponentially with rate 9q > and if q > 0, it is known that ((15)1 
(that is, SULE condition) implies the existence of a constant C^(m, I, uS) so 
that 

sup m| ,j u {m, t) < Cfy(m, P — a.s., 

i ' ' 

i.e., H\,(m, c) is dynamically localized on / (see Section 2 in |16j). 

To prove (ii) and (|15j) . it is sufficient to show strict positivity of the 
Lyapunov exponent, because in this case Lemmas 01 and hold. By using 
the multiscale analysis |28j together with (PI) and (P2), one can then follow 
the proof of Theorem 3.1 in JS] (properly adapted to IL\,(m,c)) to obtain 
(ii) and (|T5|) (details will be omitted). 
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Now the proof of (i). It follows from Furstenberg and Kesten Theorem [S] 
that, P-a.s. the Lyapunov exponent 7 m exists and is independent of u. 

Consider first the energies E ^ ±V and it will be proven that "f m (E ^ 
±V) > for all m > and for all E G a{T) w (m,c)). Let Q m (E) be as in 
the Lemma |H Put a = E - V, j3 = E + V and rename T^(E) = , 
T~ v (E) = Tm . In the present case and (3 ^ 0. 

Since the problem is symmetric in a and /3, the proof is reduced to the 
study of three cases: 

a) and T$ are elliptic ((trace T^ ] \ < 2, (trace T%p\ < 2); 

b) Tm ^ is parabolic ((trace Tm \ = 2); 

c) Tin ] is hyperbolic ((trace | > 2). 

Note that in cases b) and c) the group Q m (E) is not compact. 

Case a) Since T m ^ and Tfj^ are both elliptic, then \o\, |/3| G {rnc 2 , cy/A + tti^c 2 ) 
In this case such matrices do not commute. Since the operator 

rp{a)rp(J3) (rp{a)\-l(rn(l3)\-l 

m m \ m I \ m I 

built from two noncommuting elliptic elements is hyperbolic, it follows that 
Q m {E) is not compact. Moreover, note that 

trace {T^f = ^ - Urn 2 + A) a 2 + m 2 c 2 (4 + m 2 c 2 ) + 2 

(analogous for T^f "*) . Hence, if a 2 / 2c 2 +m 2 c 4 or /3 2 7^ 2c 2 +m 2 c 4 , then T^f ^ 
and (Tm^ 2 or and (Tm^) 2 are elliptic. Since elliptic elements have no 
fixed points in P(IR 2 ), it follows that for any x G P(IR 2 ), Q m (E) -x contains 
at least the three elements x, Tm^ • x, (Tm^) 2 ■ x or x, • x, fc^f ) 2 • x. 
Therefore, by Lemma|21 j m (E) > 0. If, on the other hand, a 2 = 2c 2 + m 2 c 4 
and (3 2 = 2c 2 + m 2 c 4 , then E = and V = cy2 + m 2 c 2 , which is one of the 
excluded pairs described in Theorem |IJ 



Case b) Suppose Tm is parabolic, that is, \a\ = mc 2 or \a\ = c\/4 + m 2 c 2 . 
First the possibility a = mc 2 will be discussed (the case a = —mc 2 is 
similar). In this case 

rw=(j 2 r), and so wr=(j 2n r). 

Denote by {ei, 62} the canonical basis of IR 2 . By taking a vector x = x\e\ + 
0:262 j and setting x for its direction, one concludes that limj^oo (Tm) n -x = 
ei . If f is a probability measure that is invariant under the action of Q m (E) , 
and if / G C^°(P(IR 2 )), by Lebesgue's dominated convergence Theorem one 
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has 

/(Si) = Jim | / ((TW) n • z) dv(x) . 

This means that f = <5gi • But the matrix T^ does not leave invariant the 
direction e\ since 

m f*n ( m 2 c 4 -(5 2 \ -13 + mc 2 , „ . 2 

Tif > ei = 1 + ei + — e 2 and ^mc 2 . 



m 2 c 2 



Thus it is proven that there is no invariant measure under the action of 
Gm{E)- Therefore, by Lemma [5] one gets "f m (E) > 0. 

Consider now the possibility a = cy/4 + m 2 c 2 (the case a = —cy/4 + 
is similar). In this case an eigenvector of 

-3 mc + \/4 + m 2 c 2 

mc — v4 + m 2 c 2 

, / mc + \/4 + m 2 c 2 A / — mc + \/4 + m 2 c 2 
is given by ui = I , 1 J . Picking v 2 = I , -1 

a vector orthogonal to «i, the matrix in the basis {^1,^2} is 



— 1 —4 — m 2 c 2 + mc\/ 4 + m 2 c 2 

-1 

Repeating the previous calculation for this case, one obtains v = 5^ . But 
Tm^ does not leave invariant the direction v\ except for (3 = or (3 = 
csj 4 + m 2 c 2 = q, which are excluded since the first condition yields E = —V 
and the second one V = 0. Thus it is proven that there is no invariant 
measure and, by Lemma |2 j m (E) > 0. 



Case c) Suppose now that Tm is hyperbolic (so \a\ < mc 2 or \a\ > 
cV 4 + m 2 c 2 ). It is sufficient to study the orbit of the eigendirections of 
Tm , namely 

(a 2 — m 2 c 4 + e^J (a 2 — m 2 c A ){a 2 — m 2 c 4 — 4c 2 ) 
2c(a - mc 2 ) 

If T^f^ is hyperbolic then the orbit of e m is infinite. Hence j m (E) > by 
Lemma|21 If T^ is parabolic, it is again case b). Finally, suppose that T^ 
is elliptic. If Tm e m 7^ e~ e , then Tm e m cannot belong to the eigendirections 
of Tm^ and its orbit is infinite. Hence "f m (E) > by Lemma EJ If Tm& m = 
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e TO e , then simple calculations lead to the equations 



1 + m 2 c 2 



B 2 

— ) (a 2 — m 2 c A + eu) = 4(m 2 c 4 — Ba) + (a 2 — m 2 c 4 — < a ) 



with e = ±1 and u = \J (a 2 — m 2 c' l )(a 2 — m 2 c 4 — Ac 2 ) ^ 0. It implies 
B 2 = 2c 2 + m 2 c 4 and a = 8 ± c\/2, which means V = c/y/2 and E = 
±c\/2 + m 2 c 2 — c/y/2. The symmetric case where one assumes that T^P 
is hyperbolic leads naturally to a 2 = 2c 2 + m 2 c 4 and 3 = a ± c\/2, which 
means V = c/\/2 and £ = ±o/2 + m 2 c 2 + c/\/2. Those are excluded pairs 
that will be discussed in the proof of Theorem 01 



Consider now the energy E = V (the case E = — V is analogous). Note 
that a = and B = 2V. First the case m > will be discussed. The two 
possible transfer matrices are 



^^(^^ 7) and T(f) 



/ m 2 c 4 -4y 2 mc 2 + 2V \ 
1 ~t~ o 



mc 



2V 



Observe that an d do not commute, and that is hyperbolic. 
It is enough to study this case for B = cV 4 + m 2 c 2 (T^f ^ is parabolic) . The 
eigendirections of T^ff 1 are 



' mc + 5V 4 + m 2 c 2 \ 
2 

V 1 / 



±1. 



The matrices and T^f in the basis {e^,^ 1 } are given, respectively, 
by 

-1 4 + m 2 c 2 



Ai 
A_i 



and 







mc 



-1 



v4 + r))-r- 



2„2 



with 
AiA_! 



1 + 



m 2 c 2 



+ 



mcy/A + m 2 c 2 



1 + 



m 2 c 2 



mc^J '4 + m 2 c 2 



Suppose that T^f ^ occurs with probability < p < 1 and IX occurs with 
probability 1 — p. Denote by n a (resp. rip) the number of times that Tffl 
(resp. Tm^) occurs in the product T^(E;n,l). Supposing, without loss of 
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generality, that Tm '(E) = T£> , one has 

A™ Q C n P(Ai,A_i) 



T£(E;n,l) 







A^ 



P-a.s., where C n is a constant and P(Ai, A_i) is a polynomial in Ai and A_i. 
Thus, 



T^(E;n,l) 
and therefore P-a.s. 



> 



A™° 




A™", Ai > 1, 



1 TL 

lm (E = V)= lim — In \\T"(E;n, 1)|| > (In A x ) lim ^ = (In Ai)p > 0. 

n^oo fl n— >oo ?i 

Now the case m = will be treated. In this case 

1 - AV 2 /c 2 2V/c ' 



T^ a) =Id 2 and T (/3) 



-2V/c 



1 



One then finds 



if V € (0, c] and 



lim \\(T^) 



(/3)\n,,l/n 



lim (T f 



(/3)Nn,|l/n 







> 1 



if V > c. Hence, if V 6 (0, c], V~ ^ c/>/2, 



7o(^ = F)= lim ^ In ||(T W ) 



(/ 3 )\ n /3|| 1 / n /3 



(1 - p) In 1 = 0, 



and 7o(£ = V) > if V > c. 



□ 



Proof, f Theorem l4|l 

By analyzing the proof of Theorems [21 and [3] observe that for V = c/y2 
(resp. V = cv2 + m 2 c 2 ) one has ^ m (Ey ^ ±cv2 + m 2 c 2 ± c/y/2) > 
(resp. 7 m (£V 7^ 0) > 0) and then the same conclusions of Theorem [21 
(resp. Theorem 0J) hold. It remains to show that 7 m vanishes at the pairs 
(V = cV2 + m 2 c 2 , £y = 0) and (V = c/v 7 ^, £y = ±cy/2 + m 2 c 2 ±c/y/2 ) . 
Note that in all cases Ey G <7(H\,(m, c)) P-a.s. . 

First it will be treated the case (V = c/y/2, Ey = —cV 2 + m 2 c 2 — cj\pl ) 
(the others three excluded cases with V = c/y/2 are similar). In this case 
one has (3 = —cs/l + m 2 c 2 and a = j3 — c\[2. The eigenvectors of Tm are 
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given by 



/ 2c - yg/3 + ey^c 2 + 2m 2 c 4 - 2\Z2c/3 \ 



V 



1 



±1, 



/ 



and by looking at the matrices in the basis given by these two vectors, the 
study is reduced to products of matrices of the following two types: 









A_ 



with A + A- 



and 



1 and k + k 



and 
1, where 











1 | V2f3 ± V4c 2 + 2m 2 c 4 - 2y/2c/3 



K± 



Moreover, 







+ mc 



2c - V2f3 ± y^c 2 + 2m 2 c 4 - 2y/2c(3 



+ 1. 



WW 



mc + (3/c 



-Id 2 



mc 



(3/c 



Therefore the proof that j m (Ey = — cy/2 + m 2 c 2 — c/y/2, V = c/y/2 ) = 
is analogous to the Schrodinger case (see the proof of Theorem 2.4 in |10| ^ . 

0). 



Now consider the excluded case (V 



cy/2 + m 2 c 2 , 



Ey 



In this 



0), then a 



Id 2 and T^T^ 



case a" = (3 2 = 2c 2 + m 2 c 4 . Since a / (3 (otherwise V 
-(3 = ±cV2 + m 2 c 2 . Noting that {T^f = {T^ ] f -- 
is hyperbolic, the proof that j m (Ev = 0, V = cV2 + m 2 c 2 ) = is again 
similar to the corresponding Schrodinger case (see the proof of Theorem 2.4 

in ma). □ 



Proof. (Theorem 15]) 

(i) It is sufficient to prove the Theorem for VP = <5q . Given A > 0, there 
exists b' > and by Lemma El there are b > and C < oo such that, 
by applying Lemma El for N = [b'T] , together with the relation 1)1 ljl for 
m = and ( = i/T, one concludes that there exists a set r2jv(A) C with 
P(Ojv(A)) < Ce~ 67vA and 

(16) \\T£(E + i/T;n,l)\\ < C 

for all to e Q\n N (X), I < n < N and E £ I v = [E v - N~ x ~ l / 2 ,E V + 
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Supposing that 

\G+JE + i/T; 1)| 2 + \G^(E + i/T; 0)| 2 > B^u) > 0, 
it follows from and 

WT^E + i/T^n^W = \\T£(E + i/T;n,l)\\ 

that 
(17) 

ma X { | GUE + i/T;n)\MG,jE + z/T;n-l)\ 2 }> B ' ! ^ 1 



2||T-(£ + i/T;n,l)|| 2 " 



Thus, replacing (|16f) and (fTTj) into (fTl]) . P-a.s. one has 

0<n<[b'T] Jlv 

for some constant C q {oj) > 0. 
If, on the other hand, 

\G+ >UJ (E + i/T; 0)| 2 + \G^(E + i/T; -1)| 2 > fl 2 ( w ) > 0, 

then one gets this estimate in the same way, but based on (|T2|) instead 
of (|13|) . Since A > is arbitrary, this finishes the proof. 

(ii) It follows from the above arguments by using Lemma |H1 □ 



Proof. (Theorem E3) 

The arguments will be a variation of [22]. Define the operator 



p := i [TD(m,c),X] = ci ( d ° 1 d Q 1 j 

with [•, •] denoting the commutator. Note that p is self-adjoint and bounded. 
Set 



X(t) = ei^(rn,c)t Xe -iJD(m,c)t &nd ^ = e ».ID( m , C ) ipe -iID(m, C )^ 



so that 



-|x(t) = i [lD(m, c),X(t)] = e ilD ( m ' c )*i [lD(m, c), X] e -^( m . c )* = p ( t ). 



Hence 



X(t) = X+ [ p(s)ds. 
Jo 



Using this relation, the boundedness ||p(i)|| = ||p|| < oo for all t, Cauchy- 
Schwarz inequality, and keeping only the dominant terms for large t, it 
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follows that for t > 1 and q S IN, there exists C q (V, m, c) > so that 



M^(m,t) = {V,\X(t)\«*) 
< C q (V,rn,c) 



<*,p(si)---p(s g )*) d Sl ---ds q 



< C q (V,m,c)\\p\\n« = K q (V,m,c)ti. 



□ 



6. Dynamical Comparison 



The aim of this section is to compare the dynamical moments M^\m, t), 
as in Definition^ for different masses and general potentials, in particular 
for the massless and massive Bernoulli cases. 



Theorem 7. Let ID(m,c) and TD(m',c) be Dime operators on Z 2 (2,(D 2 ) 
defined as in with the same potential V , and let ^> be with only one 
nonzero component. Given T > 0, there exists a constant B q > so that 



(18) 



sup 

0<t<T 



M^\m,t)-M^(m',t) 



(?), 



< B q \m -m'\c 2 T q+2 . 



Proof. Observe that for m = m' the result is immediate. Suppose m ^ m! . 
For the proof it will be assumed that m > 0, m' = and ^ = 5q (the case 
m > 0, m' > and as in the hypotheses is similar). 
For fixed a > consider the Banach space 



B n 



$ G Z 2 (2,(D 2 ) : ||$ 



i„ -supe~'^|(S+ $)| + |<<5fc,$)| <^ 



Since ID(m, c) is a bounded operator on B a , it follows that 
(19)((Ce" iID(m ' c) *^)| + \(S-,e- mim ' c)t 5+)\ < ||e~ iID(m ' c) '5+|| a e~ na 

<• — na+t||K)(m,c)||a 

For A; G IN denote by X k the restriction of the position operator X to 

the set |7T,| < k} and by Mj+ (m, i) the corresponding dynamical 

°o 



moments. Then, for all times t < 



a k 



AI 



2||ID(m,c) 
V(m,t)-M$' k (m,t) 



using (fl9|) one has 



(«), 



E 

|n|>fc 



(20) 



<Ci(g)A; 9 e- /ca+2 * l|ID(m ' c)l1 " < C^k" 
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Furthermore, it follows by DuHamel's formula that 



M^' k (m,t)-M'l>' K (0,t) 



(«).*/ 



-i J* (s£, jV{m,c)t\ x k\q e -iTD{m,c)(t-s) ^ _ ^ c jj g -iID(0,c) S( ^ ds 

+t y ((5(1", e iID(m ' c)( '- s) (lD(m, c) - 1D(0, c))e iID(0 ' c)s |X fc | <7 e- iID(0 ' c) ^^ da. 



Hence, for t < 



a k 



2||E>(m,c) 



-, using (|19|1. the fact of the operator e 



iID(m,c)t 



on £ (Z; (D ) be unitary and Cauchy-Schwarz, it is found that 



(21) M^' k (m,t)-MS' k (0,t) < C 2 (q)mc 2 k ( i +1 t e~ ka+t ^ m ^ 



r(<?): fc 



< C 2 (q)mc 2 



o 



k* +2 . 



Thus, by (JUJ) and (j2B), 



M$>(m,t)-M™(p,t) 



< B a mc 2 — n 
~ q 2\\TD(m,c 



2\\JD(m,c)\\ a 

k"+ 2 , 



a 



for all times i < 



2||D(m,c)|| Q 

Now, for each T > choose to be the smallest integer such that 



k > 



2||ID(m,c) 



a 



Therefore, for all t < T, 



M$(m,t)-M$(0,t) 



T. 



< B q mc 2 T q+2 . 



□ 



With respect to the Bernoulli-Dirac model @, the relation with m > 
and m! = gives an estimate of how, for small times and/or sufficiently 
small mass, the dynamics of the localized regime follows the delocalized one. 
In terms of the average dynamical moments A$ (m,T) defined in (|1U|) one 
has 



Corollary 1. Let (IL\,(m, c)) u&n 6e as m (0), m > and V G (0, c],V / 
c/v2j a^rf /ei VP 6e iwzi/i emZy one nonzero component. Then, for each q > 0, 



P-a.s. there is C q ^ > so that 



A (l) 



4iM 



<C^mc 2 t 5 / 2 , 



t > 0. 



Notice that the power exponent on the r.h.s. of this expression does not 
depend on q. 
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Proof. By Theorem with m > and m! = 0, it follows that for all t > 

A^Jm, t) - 4 9) J0, t)\ < T(q + 3)(V mc 2 ^+ 2 , 

with r the usual gamma function. By Theorem|5^i) there is < B„{oj) < oo 
such that 

A^,t)>B q {^~ 1 ' 2 P-a.s., 
and the result follows with C q ^ = T(q + 2>)C q ^/ B q {uj). □ 

Remark. By using Theorem one gets (for q £ IN) 

il4 9) (m,t) - M^m'.t) < k q (u,m,m',c) t q , 

but with no expression for the constant K q {uj,m,m\c). The price paid for 
the explicit dependence on the masses and light speed c in Theorem is the 
larger exponent q + 2 instead of just q. In the same way, the exponent 5/2 
in Corollary Q could be replaced by 3/2, but with no precise dependence of 
the resulting multiplicative constant on m and c. 
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